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ABSTRACT 

The small-scale magnetic field is ubiquitous at the solar surface—even at high latitudes. Erom observations 
we know that this field is uncorrelated (or perhaps even weakly anticorrelated) with the global sunspot cycle. 
Our aim is to explore the origin, and particularly the cycle dependence of such a phenomenon using three- 
dimensional dynamo simulations. We adopt a simple model of a turbulent dynamo in a shearing box driven 
by helically forced turbulence. Depending on the dynamo parameters, large-scale (global) and small-scale 
(local) dynamos can be excited independently in this model. Based on simulations in different parameter 
regimes, we find that, when only the large-scale dynamo is operating in the system, the small-scale magnetic 
field generated through shredding and tangling of the large-scale magnetic field is positively correlated with 
the global magnetic cycle. However, when both dynamos are operating, the small-scale field is produced from 
both the small-scale dynamo and the tangling of the large-scale field. In this situation, when the large-scale 
field is weaker than the equipartition value of the turbulence, the small-scale field is almost uncorrelated with 
the large-scale magnetic cycle. On the other hand, when the large-scale field is stronger than the equipartition 
value, we observe an anticorrelation between the small-scale field and the large-scale magnetic cycle. This 
anticorrelation can be interpreted as a suppression of the small-scale dynamo. Based on our studies we conclude 
that the observed small-scale magnetic field in the Sun is generated by the combined mechanisms of small-scale 
dynamo and tangling of the large-scale field. 

Subject headings: Magnetohydrodynamics - turbulence - dynamos - Sun; magnetic fields - sunspots 


1. INTRODUCTION 

Sunspots and the solar cycle are associated with the large- 
scale magnetic field of the Sun. This magnetic field is be¬ 
lieved to be produced by a large-scale dynamo, globally oper¬ 
ating in the solar convection zone. The a effect (originating 
either from helical flow and/or the Babcock-Leighton process) 
and the fl effect (shearing of the magnetic field by differential 
rotation) are bel ieved to be the cause of this dynamo in the Sun 
(see reviews bv lCharbonneaul2014tlKarak et al.ll2014ll . How¬ 
ever, there is another type of magnetic field that exists at small 
scales and even in the quiet phase of the Sun. This field is in¬ 
termittent and fluctuating and exists in mixed polarity even at 
the resolution limits of present day instruments. This field 
is known as small-scale, fluctuating, turbulent, or internet¬ 
work field—as it exists in the internetwork regions. Although 
the e xistence of this field has been known for a long time 
(e.g.. [Howard & StenflolflOT^ iFrazier & Stenflol 11972 ). the 
detailed understanding of its nature and origin has remained 
elusive. The main reason for this is that the spatial charac¬ 
teristic scale of this field is shorter than the resolution limits 
of present day instruments and the weak polarization signal 
makes the measurement of this quiet-Sun magnetic field dif¬ 
ficult (e.g.. lde Wiin et al1l2009h . Interestingly, even the inter¬ 
network field—the weakest component of solar magnetism— 
having an unsigned flux of ~ 10®® — 10®® Mx, contributes 
about 10^® Mx of magnetic flux to the solar surface per day, 
which is at least four orders of magnitude more than that 
from the bipolar active regions (which contain an unsigned 
flux of ~ 5 X 10^° — 5 X 10^^ Mx) during solar maximum 
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(ISchrijver & Harve^ll994IZhou et al.ll2013h . In comparison, 
both the network regions with a flux of Ri 10®®^ — 10®® Mx 
and ephemeral regions with 10®® — 5 x 10®® Mx contribute 
rough W two orders of magnhude more flux th an active re¬ 
gions ^agenaareT^lIlQS |GolIl]iLiDl20l3)- Therefore, 
the small-scale magnetic field is potentially important in un¬ 
derstanding many unsolved problems in solar physics, parti c- 
ularly s olar atmospheric h eating (e.g.. lSchriiver et ^119981) . 

Since iBatchelon (119501) it is known that a sufficiently ran¬ 
dom three-dimensional velocity field can amplify a seed 
magnetic field via repeated random stretching, bending and 
folding—without requiring net helicity. This is the essence 
of the small-scale or local dynamo. The correct mathe- 
ma tical description o f such dy namos was develop e d late r 
by iKaza ntsevI (11 96^ : see also iChildress & Gilbert (119951) . 
iPetrovav & Szakalvl (119931) were the first to realize the ne¬ 
cessity of a small-scale dynamo t o explain the magne tic flux 
budget in the photosphere (also see lDurnev et alJl993l) . How¬ 
ever, one may put forward alternative arguments in which 
this quiet-Sun magnetic field is the result of a global large- 
scale dynamo; the shredding of large-scale magnetic field 
can produce small-scale magnetic field by cascading mag¬ 
netic energy to small-scales and the deca y of active regions 
may give rise to small-scale magn etic field (ISpruit et al.ll 19871 ; 
Ide Wiin et alJ2005l ; lStenflol201^ . However, none of these ar¬ 
guments may necessarily be correct because the small-scale 
magnetic field does not show solar cycle dependence, as 
it does not have significant correlation with the large-scale 
global magnet ic cycle, and it does not have much latitudinal 
variation (e.g..lHagenaar et al.ll2003l; fSanchez Almeidall2003l; 

iLites et aHl2008l;lLitesll201lHBuehler et al.ll2013l;lJin & Wand 
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l2015bl) . Furt hermore, using the MDI data set of the entire so¬ 
lar cycle 23. lJin et'dl (1201 Ih have shown that the cyclic vari¬ 
ation of both number and total flux of small-scale magnetic 
elements in the range ~ (3-30) x 10^® Mx are anticorrelated 
with the sunspot number. Subsequently. iJin & Wand (I2015al) 
have studied the internetwork magnetic field within the flux 
range of about Mx and And no correlation with the 

solar cycle. All these observational results support the small- 
scale dynamo scenari o as the source of sm all-scale held in 
the Sun. We note that lSokoloff et al.l (12015h argued in favor 
of the operation of a small-scale dynamo in the Sun to ex¬ 
plain the observational findings that the number of anti-Hale 
sunspot groups becomes maximal during solar minimum and 
that these groups do not show the latitudinal dependence as 
Hale sunspots'. 

The complex nature of the stratified convection zone with 
very low values of molecular viscosity, diffusivity, and mag¬ 
netic Prandtl number (Pm) make the theoretical modeling of 
the small-scale held challenging at present. Furthermore, 
at small Pm, the critical magnetic Reynolds number, 
needed to excite the small- scale dynamo increases at first 
linearly with deceasing Pm ( Rogac h evskii & Kleeorin|ri997l: 
iBoldvrev & Cattaneol l2004t Schekochihin et al.l 120051) until 
P£ s aturates for Pm 1 (llskakov et al.ll2007l: iBrandenburgI 
1201 ll) . which makes the small-scale dynamo at the solar 
value of Pm ^ 1 0~® possible. Beside these difficulties sev- 
eral authors (e.g.. lCattaneol[T9^ lEmonet & Cattaneoll20()n) 
tried to model the small-scale dynamo in a local Cartesian 
box. Although they found in some cases small-scale dy¬ 
namo action, these simulations are highly ideal i zed an d far 
fro m the real Sun. Subsequently, IV ogler et^ (12005 ) (see 
also iVogler & Schusslefll2007t iPietarila Graham et al.l 2010t 
lRempelll2014ft have gone one step further to make these mod¬ 
els more realistic by including radiation and ionization. How¬ 
ever, these simulations are local in the sense that they cannot 
study the effects of the global dynamo-ge nerated large-scale 
field on the small-scale field. Very recently. fHotta et ^ (12014 
120151) have studied small-scale dynamo action in the Sun us¬ 
ing high-resolution MHD simulations. Again, however, these 
simulations did not consider the large-scale dynamo. 

The purpose of the present study is to explore the effects of 
a global dynamo-generated cyclic large-scale magnetic field 
on the small-scale field. To this end, we perform a set of 
simulations of a turbulent cyclic dynamo that produces a mi¬ 
gratory magnetic field similar to that of the Sun. To capture 
the essential features of an ail dynamo, in particular its mi¬ 
gratory properties, we employ helically forced turbulence to¬ 
gether with imposed uniform shear across the domain. These 
two ingredients provide a prototype of an aH dynamo, as it is 
often modelled in mean-field theory. Earlier turbulence sim¬ 
ulations of this type have confirmed the presence of cyclic 
solutions with dynamo waves whose direction of propaga- 
tion depends on the sign and magnitude of shear a nd helicity 
(iKapvlk & Brandenburgl2009t iKarak et al.ll2015bh . 

We do not make an attempt to model the real Sun, whose 

* One may argue tliat tire anti-Hale sunspot groups could be 
formed through the negative effective magn etic pressure instability 
(Brandenburg et al.l201 lLl2013l : IKemel et al.l2013l) because, for this to work, 
we only need a moderately weak magnetic field, which could be provided 
by the small-scale dynamo operating continuously in the Sun. The anti-Hale 
grou ps would then represent just the tail of a wide distribution of tilt an¬ 
gles (McClintock et alj|2014ll . Within the rising flux tube picture, anti-Hale 
groups near the equator could alternatively be explained as tubes buffeted 
from the other hemisphere. 


underlying dynamo mechanism is still not well understood 
(lBrandenburg|l2005l: ICharbonneaull20 1 4) . It is not even com¬ 
pletely clear whether the solar dynamo is dr iven by con¬ 
vective motions or by m agnetic instabilities (IBrandenburgI 
119981: iTobias et al.l I2()l ll) . and whether meridional circula¬ 
tion plays an important role for causi ng the equatorward mi¬ 
gration of the Su n’s magnetic field (IChoudhuri et al.l [19951: 
iHazra et al.l20l4) . Furthermore, we are facing even more fun¬ 
damental challenges, for example the question of catastrophic 
quenching, the dependence of mean field strength on mag¬ 
netic Reynolds number, and the morphology of the field in 
the form of tubes on the one hand, and the question of mag¬ 
netic buoyancy, flux tube storage, and sunspot formation on 
the other. Significant progress has been made by considering 
these aspects separately, since in models with all effects in¬ 
cluded it is difficult to associate particular features to given 
changes in the model. Thus, in our model, we focus on 
the aspect of large-scale and small-scale dynamos only, i.e., 
we do not consider the question of magnetic buoyancy, flux 
tube storage, and sunspot formation. Consequently, we con¬ 
sider a model without gravity and convection and consider 
turbulence to be driven by a stochastic forcing. Moreover, 
our simple model then avoids the vastly different time and 
length scales of global and local dynamos in realistic settings. 
Capturing them within a single model r emains a challeng¬ 
ing task in global convection simulations (IVo gler et a Hl2005l: 
iKMpvla et al.ll201!^ iKarak et al.ll2015al: iHotta et al.ll2015l) . In 

the present paper, our aim is to address the basic question 
whether and how the small-scale magnetic field is correlated 
with the dynamo cycle in a simple realization of a turbulent 
dynamo that combines the physics of a small-scale dynamo 
with that of an ail dynamo. 


2. MODEL 

In this work we do not consider convection. Therefore, 
gravitational stratification and thermodynamics are neglected. 
Instead, our model consists of an isothermal, compressible 
fluid where the pressure is given by p = c^p, with Cs being 
the isothermal (constant) sound speed and p the density. We 
sustain turbulence by supplying energy to the system through 
a stochastic forcing function in the momentum equation. As 
we do not include the stratification and global rotation of the 
Sun, we make this forcing helical to produce an a effect in 
our model. Furthermore, to mimic the large-scale shear of the 
Sun and to produce an it effect, we impose a uniform shear 
flow in our model. We solve the hydromagnetic equations in 
a cubic domain of size (27r)®. To connect with spherical ge¬ 
ometry, we can imagine our simulation box being located in 
the northern hemisphere of the Sun such that the x axis cor¬ 
responds to the radially outward direction, the y axis to the 
toroidal direction, and the z axis t o the d irection of increasing 
latitude; see, e.g., iKapvla et al.l (120131) . who used the same 
correspondence in their Section 3.6. Within this framework, 
the momentum, continuity, and induction equations become 
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where D/Dt = d/dt + {U + ■ V is the advective 

time derivative, = {0,Sx,0) is the imposed uniform 

large-scale shear flow with S being the constant shear rate, f 
is the forcing function that we describe later, i/ is the constant 
kinematic viscosity, and S is the traceless rate of strain tensor 
given by Sy = ^{Uij + Uj^i) — V • U (the commas 
denote partial differentiation with respect to the coordinate j 
or i). The contribution of to the viscous stress is small 
compared with that of U and is neglected. We note that the 
induction equation is written in terms of the magnetic vector 
potential A, such that V x A = B is the magnetic field. Fur¬ 
thermore, J = X B is the current density, and rj is the 

constant microscopic diffusivity 

The forcing function f = fix A) in Equation ([T]) is 
helical and random in time. It is defined as f{x,t) = 
exp[ife(f) • x -|- i(/)(f)]}, where x is the posi¬ 
tion vector. At each timestep, the wavevector fc(f) randomly 
takes any value from all possible wavevectors in a certain 
range around a given forcing wavenumber fcf. The phase 
—TT < (j)it) < TT also changes randomly at every timestep. 
On dimensional grounds, we choose N = /oCs(|A;|cs/(5f)^/^, 
where /o is a non-dimensional forcing amplitude. The trans- 
verse helical wave s are produced via Fourier amplitudes 
(iHaugen et al.ll2004h 

= with ( 4 ) 

where cr is a measure of the helicity of the forcing; 
for positive maximum helicity, cr = 1, = 

(fc X e) / yjk'^ — [k ■ e)2, where e is an arbitrary unit vector 
not aligned with fc. Note that |/fcp = 1 and /fc • (ifc x /fc)* = 
2crfc/(1 + cr^). 

We solve the above equations with the PENCIL CODE^ 
in—TT < x,y, z < tt with boundary conditions that are 
shearing-periodic in the x direction and periodic in the y and 
z directions. The numerical time step Af is computed as the 
Courant time step. As initial conditions we take tt = In p = 0 
and A consists of small-scale low amplitude (10“^) Gaus¬ 
sian noise. In our nondimensional unit system, we have 
Cs = Pa = yo = 1, where po = {p) is the volume-averaged 
density, which is constant in time because of mass conserva¬ 
tion. 

To interpret our results in terms of a mean-field picture, we 
define the mean-field dynamo number D = CaCn with 

C'a = ao/pTofci and Cq = \S\/r]Tokl, (5) 

where pxo = Vw+il is the theoretically expected total (turbu¬ 
lent plus microphysical) magnetic diffusivity, pto = 
ao = —t{u} ■ u) /3, T = (urmsfcf)~^, aud kf is the mean 
forcing wavenumber (iBlackman & Brandenburg|l2002h . The 
other usual diagnostic parameters, i?m, Re (fluid Reynolds 
number), and are defined as 

Re = Bra — ^/V-j (A) 

where Urms = is the rms value of the velocity in 

the statistically stationary state with (•) denoting the average 
over the whole domain and over time, and i?eq = Wrms is 
the volume-averaged equipartition held (we recall that in our 
units, poip) = 1). Furthermore, the large-scale (mean) field 

^ http://github.com/pencil-code 
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Figure 1. Different regimes in magnetic Reynolds and dynamo numbers 
space to be studied in the present work. 

will be defined as a horizontal average and therefore the small- 
scale held will be defined as the residual of the total and the 
large-scale quantities viz, 

lP=W-B^. (7) 

Here again overbars denote xy-averages. 

3. RESULTS 

We perform a set of simulations in four different regimes 
in the parameter space D). We define the Regimes I, II, 
III, and IV as marked in Figure [T] 

3.1. Small i?m ond small D: Regime I 

In this regime, as both the Reynolds and the dynamo num¬ 
bers are small, neither the small-scale dynamo nor the large- 
scale dynamo are possible. Therefore, any seed magnetic held 
decays with time. In Table [T] we show a typical run; R1 at 
V = 0.01, p = 0.008, (7 = 1 (fully helical forcing of posi¬ 
tive helicity), and S = —0.05, so = 0.55 and D = 4.1, 
resulting in a slow decay at a rate O.OSurmsfci- 

3.2. Small and large D: Regime II 

Next, we consider a set of parameters in Regime II. 
Run R2a is a typical simulation in this regime which is sum¬ 
marized in Table [T] In the two top panels of Figure |2] we 
show the X- and y-components of the mean held, and 
By. We see that the large-scale field is oscillatory with dy¬ 
namo waves propagating in the positive z direction. This is 
indeed expected because we force the flow with positive he¬ 
licity, which results in a negative a effect, and since shear 
is also negative, the dynamo number is positive, resulting in 
dyn amo wave propagation in the positi ve z direction (see, 
e.g.. [Brandenburg & Subramanianll2005[) . The toroidal held 
(By, streamwise direction) is about 10 times stronger than the 
poloidal held (Bx, cross-stream direction); see Figured This 
is because of the fact that the H effect, as measured by Cq, is 
much larger than the a effect, as measured by Cq,; see Table[T] 

The next panel of Figure |2] shows the energy density of the 
_ 2 

large-scale held, B , and the last panel that of the small-scale 
field, b^, as computed from Equation (|7]). An interesting fea¬ 
ture to note is that the small-scale held is spatially and tempo¬ 
rally correlated with the large-scale one. We also see that is 
much smaller than B . The time series of B^ and shown 
in Figure [3 show a tight correlation. Moreover, the small- 
scale held becomes significant only when the large-scale field 
does. In the early phase, this small-scale field grows linearly 
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Table 1 

Summary of the runs. 


Input Output 


Regime 

Run 

grid 

ki 

V 

Pm 

|S| 

Urms / Cs 

Pm 

Px 

By 

B 

b 

Cc 

Cq 

D Covr{By,b^) Dynamo? 

I 

R! 

48^ 

5.1 

10"^ 

1.25 

0.05 

0.023 

0.55 

0.000 

0.00 

0.00 

0.000 

0.78 

5.3 

4.1 

— 

No dynamo 

II 

R2a 

48^ 

5.1 

10"^ 

1.25 

0.1 

0.022 

0.55 

0.165 

1.40 

1.41 

0.274 

0.75 

10.6 

7.9 

0.77 

LSD 

II 

R2b 

96® 

3.1 

5 X 10-4 

0.42 

0.05 

0.066 

17.7 

0.246 

1.34 

1.38 

0.662 

1.38 

6.1 

8.4 

0.42 

LSD 

II 

R2c 

96® 

3.1 

5 X 10-4 

0.50 

0.04 

0.063 

20.1 

0.298 

1.50 

1.54 

0.711 

1.58 

5.2 

8.2 

0.48 

LSD 

II 

R2d 

48^ 

5.1 

10-2 

1.25 

0.15 

0.022 

0.54 

0.202 

2.57 

2.61 

0.418 

0.71 

15.9 

11.4 

0.91 

LSD 

II 

R2e 

00 

5.1 

10-2 

1.25 

0.2 

0.039 

0.97 

0.294 

4.56 

4.57 

0.671 

1.14 

18.9 

21.6 

0.76 

LSD 

II 

R2f 

48^ 

5.1 

10-2 

1.25 

0.25 

0.039 

0.97 

0.271 

5.17 

5.18 

0.699 

1.11 

23.6 

26.2 

0.70 

LSD 

II 

R2g 

00 

5.1 

10-2 

5.00 

0.08 

0.021 

2.06 

0.308 

7.17 

7.18 

0.862 

2.00 

23.7 

47.4 

0.62 

LSD 

II 

R2h 

00 

5.1 

10-2 

5.00 

0.1 

0.021 

2.06 

0.282 

8.06 

8.07 

0.884 

1.98 

29.7 

58.8 

0.50 

LSD 

III 

R3 

144^ 

3.1 

10-^ 

5.0 

0.0 

0.054 

86.1 

0.076 

0.11 

0.15 

0.526 

0.00 

0.0 

0.0 

— 

SSD 

III 

R3' 

288® 

3.1 

10-3 

5.0 

0.0 

0.054 

86.1 

0.076 

0.11 

0.15 

0.528 

0.00 

0.0 

0.0 

— 

SSD 

IV 

R4a 

144^ 

3.1 

10-3 

5.0 

0.05 

0.057 

90.7 

0.216 

1.83 

1.88 

0.882 

2.07 

8.0 

16.6 

-0.57 

SSD -b LSD 

IV 

R4a' 

288^ 

3.1 

10-3 

5.0 

0.05 

0.057 

90.6 

0.235 

1.88 

1.90 

0.937 

2.14 

8.0 

17.2 

-0.56 

SSD -b LSD 

IV 

R4b 

1443 

3.1 

10-3 

5.0 

0.1 

0.065 

104 

0.155 

1.42 

1.45 

0.835 

1.41 

14.1 

19.8 

-0.53 

SSD -b LSD 

IV 

R4c 

1443 

3.1 

4 X 10-4 

5.0 

0.05 

0.059 

234 

0.226 

1.72 

1.74 

0.983 

2.38 

7.9 

18.8 

-0.69 

SSD -b LSD 

IV 

R4d 

1443 

3.1 

2 X 10-4 

1.0 

0.05 

0.069 

110 

0.231 

1.65 

1.68 

0.831 

1.87 

6.7 

12.5 

-0.65 

SSD -b LSD 

IV 

R4e 

1443 

3.1 

10-4 

0.5 

0.05 

0.069 

110 

0.230 

1.56 

1.58 

0.828 

1.92 

6.6 

12.7 

-0.63 

SSD -b LSD 

IV 

R4f 

1443 

3.1 

4 X 10-4 

5.0 

0.05 

0.056 

225 

0.211 

1.44 

1.47 

1.005 

1.99 

8.23 

16.4 

-0.62 

SSD -b LSD 

IV 

R4g 

1443 

3.1 

10-3 

5.0 

0.3 

0.170 

271 

0.059 

0.53 

0.53 

0.486 

0.14 

16.4 

2.3 

0.01 

SSD -b LSD 

IV 

R4h 

1443 

3.1 

10-3 

5.0 

0.2 

0.162 

258 

0.062 

0.39 

0.40 

0.456 

0.16 

11.5 

1.9 

0.14 

SSD -b LSD 


Note. — The /cf is in unit of ki = I'k = 1. Except for Run R3 for which cr = 0, for all runs we have a = 1. For Runs R2e and R2f we have /q = 0.02 
while for all other runs /o = 0.01. The B is the temporal mean in the statistically stationary state of the large-scale magnetic field over the whole domain, 
Brnis = + {By)l + (^ 2 )y>x{£’ normalized by Beq. Similarly, Bi = )^/^/Beq for i = x and y, and b = Here, Con'(By, b‘^) 

is the linear correlation coefficient r between Bf. and b^. SSD and LSD stand for small-scale and large-scale dynamos. Bold fonts in the left column indicate 


regimes of our primary interest. 



Figure!. Results from Run R2a in Regime II: Butterfly diagrams of Bx, 
By,B, and as functions of 2 and t, normalized by the diffusive time 
scale Tdiff = {ki being the smallest possible wavenumber in the 

box). The field is always normalized by Beq. The initial phase for £/rdiff < 
30 is not shown, because the magnetic field is still too weak. 

in time, contrary to the exponential growth when the field is 
generated through the small-scale dynamo. 



Figure 3. Results from Run R2a: Time series of B^/B^q (solid black 
line), b^/Bgq (red dashed), and B^^{t)/B^^ (blue dash-dotted). Note that 
is rescaled by a factor of 15. 


A scatter plot between By and for all ^ values shows a 

clear positive correlation with ~ qBy-, see Figure |4] It is 
interesting to compare this with results from the second or¬ 
der correlation approximation, which yields q ~ when 
Rm 1 and the field i s sufficiently weak; see Sect. 3.11 of 
iKrause & Radleii (119801) . However, here we find q ^ 0.03, 
which is about one order of magnitude smaller than 
(~ 0.3). One reason may well be that here the large-scale 
field is not imposed, as assumed in the analytic theory, but is 
the result of the large-scale dynamo. 

Our results demonstrate that in the present regime of small 
Rm, the small-scale field is the result of tangling of the large- 
scale field and it is not due to the small-scale dynamo. We 
recall that the small-scale field developed in the global con- 
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Figure 4. Results from R2a: Scatter plot between and 

at f/rjiff > 40 (nonlinear phase) for all z. The inset displays the same 
variation, but in linear scale. The linear correlation coefficient r between 
these two energies is 0.77. 
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Figure 5. Same as Figure|4] but from Run R2c. The red dashed line shows 
Equation —the theoretical variation. 


vection simulation of iRacine et al.l (1201 ll) in spherical geome¬ 
try also shows a positive correlation with the cyclic large-scale 
magnetic field. This suggests that the small-scale field in their 
simulation is also the result of tangling of the large-scale field 
rather than the local dynamo. 

We have performed additional simulations in this regime 
by changing the parameters (/cf, v, rj, or 5); see Runs R2b- 
R2h in Table [T] We notice that in Runs R2a, R2b and R2c, 

the large-scale magnetic field (B) has not changed much as 
D is similar in these runs, but and the small-scale field 

(b) have changed largely. The increase of i?m in these runs 
increases the tangling. Therefore, the small-scale magnetic 
field increases with increasing R^. Another interesting re¬ 
sult we notice is that, unlike in the previous case at small i?i„ 
(Figure |4]i, tangling does not increase all the way with the 
large-scale field. To demonstrate this we show a scatter plot 


between and for Run R2c—the case with the strongest 

7?in in Figure |5] Here we clearly observe the saturation of 

—2 

small-scale field. We also notice that when By just exceeds 
there is a small decrease of 6^—this is caused by the fact 
that, at this field strength, the flow is significantly suppressed 
which reduces the tangling. Then for By > B^^, remains 
roughly constant. 

The magnetic fluctuations produced in our simula- 
tions can be compared with the analytical theory of 
iRogachevskii & KleeorinI (I2007h who have studied the mag¬ 
netic fluctuations generated by tangling of the mean magnetic 
field through velocity fluctuations. Although in our simula¬ 
tions the values of i?m are rather small, we may try to com¬ 
pare with their analytical expression; 

^ + Cl /(c2/3), (8) 


with coefficients ci = 1/12, C 2 = 4, and /3 = By/B^q, 
as well as the quenching function f{f3) = 6 — 3Ai°^(/3) — 


^ 2 °^ (/3) where (/3) and A'^’ (/3) are given by 


l(0) 


( 0 ), 




tan ^ 


/33 


—/3^ lni?ni — 2/?^ In 


(3 + 5/?^)- 

l + /3^ 

i + /32Vic; 


). 




2 tan ^ /3 
5 f 


(9 -h 5/3^) -h ^ 


-/3^ In 


2/3^ In 


1 + ^^ A 

l + /32v^; 


(9) 


The above Equation ([8ll is obtained by ignoring the effects 
of convection, which corr esponds to putting the pa r amete r 
a* = 0 in Eq. (A21) of [Rogachevskii & KleeorinI (I2007h . 
The superscripts (0) indicate that these values are not affected 
by the mean magnetic field, so they are the kinematic values, 

hence Beq^^ is the unquenched value of and 6^ is the 
contribution from the small-scale dynamo, which is equal to 
zero as in this case there is no small-scale dynamo. We com¬ 
pute 6^ from Equation ([^ and display it as a dashed line in 
Eigure|5] We see that the analytical theory, which is valid for 
i?ni ^ 1 and Re ^ 1, is not far from our simulation data, 
which are at much smaller i?„i (« 20) and Re (Ri 40). Our 
result generally confirms the theoretical prediction of the sat¬ 
uration of small-scale magnetic field at By > B^q. 

Returning to the other runs in this regime (R2d-R2h), we 
notice that is not much different in these runs, but the 
values of D have increased systematically from Run R2d 

to Run R2h. The small-scale field (6) has also increased 
systematically in these runs. This clearly demonstrates that 
the tangling increases not only with i?ni, but also with the 
large-scale field itself. However, in all these runs, except in 
Run R2d, we again observe the saturation of small-scale field 
for By > Beq. 


3.3. Large Rm and vanishing D: Regime III 

We perform a simulation with <7 = 0 (non-helical), and 
5 = 0. Run R3 in Table [T] represents this simulation. When 
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Figure 6. Results from Run R3 in Regime III: Butterfly diagram of the en- 

_2 

ergy density of the large-scale magnetic field B (top) and time series of the 
small-scale field 6^ at 2 = 0 (bottom). 

Ca = 0 and Cn = 0, we should not expect a large-scale mag¬ 
netic field to develop. We have chosen such setup to exclude 
the effect of the large-scale magnetic field in this simulation. 
The top panel of Figure |6] shows the energy of the large-scale 
_2 

field B . We do not see any large-scale structures compara¬ 
ble to the system size. In any case, the small-scale dynamo is 
expected to be excited as in this case the magnetic Reynolds 
number ^ 86 is sufficiently large for large values of the 
magnetic Prandtl number. Pm = 5, to excite the small-scale 
dynamo. This is indeed confirmed by the present results. 

In the bottom panel of Figure |6] we show the time series 
of the energy density of the small-scale field. Initially this 
field grows exponentially in time and then remains constant. 
The characteristic length scale of this field is, as expected for 
a small-scale dynamo, smaller than the forcing scale of the 
turbulence, as can be seen from the magnetic energy spectrum 
in Figure|8](thick dashed curve). 

3.4. Large and large D: Regime IV 

Our first two models in Regime IV are similar to Run R3, 
but with finite shear and helicity such that the large-scale dy¬ 
namo is excited. We begin with S = —0.05 and ct = 1. In 
Table [H we refer to this run as R4a. The magnetic butter¬ 
fly diagrams of various quantities are displayed in Figure [T] 
Again, we see dynamo waves of the large-scale field com¬ 
ponents propagating in the positive z direction. The mean 
toroidal field (By) is here about 10 times stronger than the 
mean poloidal field (Bx). 

In Figure 0 we show for Run R4a' (a high resolution 
version of Run R4a) kinetic and magnetic energy spectra, 
Ek {k,t) (thin solid and dotted) and Em {k,t) (thick solid and 
dotted), during magnetic maximum and minimum, respec¬ 
tively. For comparison we also show the spectra for a small- 
scale dynamo from Run R3' (a high resolution version of 
Run R3). For both dynamos, Pif (fc, t) has a peak at the forc¬ 
ing wavenumber which is around 3 in these cases, followed 
by a small inertial range, that is slightly steeper than 
and finally the energy dissipation at small scales. Flowever, 
Pm(^, t) shows different behaviors in the two cases. For the 
LSD there is an intermediate peak slightly below the forcing 
scale (k ~ 4) and then there is a maximum at the largest 
scale (k = 1), as expected for the large-scale dynamo. During 
magnetic maximum, EM{k, t) is larger at large scales (com¬ 
pare thick solid and dotted curves), but it is slightly smaller at 
small scales (k > 20). However, during magnetic maximum 
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Figure?. Results from Run R4a in Regime IV: Butterfly diagrams of Bx, 
By,B^, andP. 



k 

Figures. Evolutions of the energy spectra EK(k,t) (thin curves) 
and EM{k,t) (thick) for Run R3^ (SSD, dashed curves) and Run R4a^ 
(SSD-hLSD, solid and dotted lines correspond to magnetic maxima and min¬ 
ima, respectively). All spectra ai'e temporally averaged excluding the initial 
growing phase. 


EK{k,t) is slightly smaller at larger scale (k < 3; compare 
thin solid and dotted curves). 

Owing to the cyclic variation of the large-scale field we see 
a weak modulation of the energy of the small-scale field, as 
seen in the lower two panels of Figure |7] In fact, 6^ shows 
a weak anticorrelation with This is evident from Fig¬ 
ure |9l where time series of By and 6^ at z = 0 are shown. 
From Figure|9]we note that hR reaches significant levels much 
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Figure 9. Results from Run R4a: Time series of By (solid black line), 
(red dashed), and (blue dash-dotted). Note that is rescaled 

by a factor of two. 
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Figure 10. Results from Run R4a: A scatter plot between B^/B^^, and 
6^/Bgq at different values of z coordinate. The red dashed line shows the 
6^ oc By dependence, whereas blue dotted line represents Equation m- 
The data points for initial growth phase (f/rdiff < 0.6) are not used for this 
plot. The linear correlation coefficient r is marked on the plot. 

faster than the large-scale energy. This implies that initially 
is primarily a result of the small-scale dynamo, which grows 
faster than the large-scale dynamo. Later, with the develop¬ 
ment of the large-scale field, increases further. Eventually, 
in this run, the saturated value of is much larger than that 
in Run R3 (which is basically the same as in this run, but only 
the small-scale dynamo is excited because of taking <7 = 0, 
and S' = 0; see Table[T]i. 

Considering only the nonlinear stage {By > Beq), we pro- 
duce in Figure [TO] a scatter plot between By and b^ for all 
grid points along the z direction. This shows an anticorre¬ 
lation obeying 6^ ; see the red dashed line. The 

linear correlation coefficient between B^ and b^ is —0.57. It 
is not straightforward to understand the origin of such vari¬ 
ations; the small-scale magnetic field developed in this sim¬ 
ulation can be read as = b^sSD + ^^tang. where b^sSD 
results purely from the small-scale dynamo, while b^tang re¬ 




sults from the tangling of the large-scale field. Unfortunately, 
separating these two quantities from a given simulation is not 
an easy task because both quantities are affected by the large- 
scale field B. In the initial phase, when the large-scale field is 
not yet fully developed, the small-scale field is mainly a con¬ 
sequence of the small-scale dynamo. Thus, in this phase we 
have b^ Ri b^ggp. Then, with increasing B^, the small-scale 
field is increasing further because of the contribution from the 
tangling of the large-scale field. However, when B^ exceeds 
Bgq, b^tang tends to saturate. (Such behavior was seen in ear¬ 
lier simulations of Section [32] showing that, when the small- 
scale dynamo was not operating, the small-scale field gener¬ 
ated due to tangling of the large-scale field increased with the 
mean-field first and then saturated; see Figure |5]) Therefore, 
in the nonlinear phase of the present simulation, the anticorre¬ 
lation results from the fact that the small-scale dynamo is sup¬ 
pressed by the large-scale magnetic field through the Lorentz 
force acting on the flow. However, when we look at the en¬ 
ergy of the turbulence, B‘l^{t) (= v?), as shown in Figure|9] 
we see no anticorrelation between the flow and the large-scale 
field. This is possibly because our turbulence is continuously 
forced by a forcing function which could break the anticorre- 

lation between By and u . 

Keeping in min d _ that the theory of 

iRogachevskii & KleeorinI (I2007^ was developed for magnetic 
fluctuations generated due to tangling and a small-scale 

dynamo represented by a constant factor b'^^ \ which is 
not the case here, we now try to fit our data by choosing 
different values of p and q in Equation (|8]l. We find that, 

when = 1.02i7eq°\ ci = 4.0, and C 2 = 0.5, there is 

_2 

moderate agreement, but only for 2 < By/Bl^ < 10; see 
blue dotted line in Figure [TOl 

We have performed a few more simulations in this regime 
(see Runs R4b-R4h in Table]!]) at different values of R^, Pm, 
or D. We observe that, as long as the large-scale field is above 
Peq. the small-scale field is anticorrelated with the large-scale 
field. This is the case for Runs R4a-R4f. However for Runs 
R4g and R4h we see that the large-scale dynamo is much 
weaker (compare the D values) and consequently the large- 
scale field is much weaker than Peq so that it cannot affect 
the small-scale field production significantly. Therefore, we 
do not see any anticorrelation between small- and large-scale 
fields in this case. In Figure [H] we show the time series of 

By and b^ and in Figure ]T2]the corresponding scatter plot for 
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Figure 12. Same as Figure[T^but results from Run R4g. The data points for 
initial growth phase <0.3) are not used for this plot. 


Run R4g. Here we do not see any anticorrelation between the 
small-scale and the large-scale fields. Another feature to note 

is that in the earlier Runs R4a-R4f the small-scale field b was 
more than 82% of B^q (see Table [T]i, whereas in Runs R4g- 
R4h the small-scale field is only about 48% of B^q- We would 
expect the opposite because the small-scale dynamo should be 
much stronger as i?i„ is larger in the latter two runs. There¬ 
fore, larger values of b in Runs R4a-R4f compared to Runs 
R4g-R4h imply that a significant contribution of small-scale 
field in Runs R4a-R4f is coming from tangling. 

On th is occasion we mention that, using the same setup 
as here, iHubbard et al.l (1201 ll) have demonstrated that, after 
some initial saturation phase, the afl dynamo mode is “frat- 
ricided” by a separate mode and eventually the initial afl 
model changes to an type dynamo, which is sustained by 
the helical turbulence alone. Some of the simulations, particu¬ 
larly at larger and large D (Runs R4a-R4f), display such 
behavior. Therefore, analyzing these simulations for longer 
durations is problematic. In this connection it maybe use¬ 
ful to note that the typical wall clock time per timestep per 
mesh point for simulations with 144^ meshpoints using 48 
processors is about 0.07 ps. Therefore, for example, to pro¬ 
duce results for Run R4a (presented in Figure|7]i it took about 
24 hours, while for R4g (Figure fTTI) it took 43 hours. As ex¬ 
amples, for Runs R2a and R4g, Tdifi /At ~ 2400 and 410000, 
respectively. 

4. SUMMARY AND CONCLUSIONS 

We have performed a set of simulations in different param¬ 
eter regimes of Reynolds and dynamo numbers to explore the 
effect of the cyclic variation of the dynamo-generated large- 
scale magnetic field on the small-scale field. Our ultimate 
motivation is to shed some light on the solar cycle depen¬ 
dence of the small-scale magnetic field at the solar surface. 
In the Sun, both global and local dynamos are expected to op¬ 
erate at the same time and can in principle occur even in the 
same location. To study this problem in detail, it is necessary 
to develop a model in which both dynamos are operating. We 
have therefore considered a simple turbulent dynamo in which 
the fluid is helically forced and large-scale linear shear is im¬ 
posed. This simple setup captures the essential mechanism 


of an afl dynamo believed to be operating in the Sun. At 
the same time, the local dynamo can be excited in this setup 
when i?„i exceeds a critical value. The main advantage of 
considering this type of simple setup, rather than considering 
more sophisticated global convection simulations in spherical 
geometry (Hotta et al., private communication) are an afford¬ 
able computational expense and that we can easily change the 
large-scale dynamo parameters to explore different regimes of 
these two dynamos. 

In Regime III, when the large-scale dynamo number is 
small but the Reynolds number is large, we observe only the 
small-scale dynamo in which the generated small-scale field 
does not have any temporal variation. This situation does not 
apply to the Sun as large-scale magnetic field is also gener¬ 
ating by the global dynamo. In Regime IV (Section 13.41) . 
when both the Reynolds and dynamo numbers are large^, 
we observe both dynamos—the large-scale dynamo produc¬ 
ing an oscillatory large-scale magnetic field and the local dy¬ 
namo producing a small-scale field. This situation is relevant 
to the Sun. In our simulations, we observe that the small- 
scale field is almost uncorrelated with the dynamo cycle as 
long as the large-scale field is smaller than the equipartition 
value of turbulence. However, when the large-scale field is 
stronger than the equipartition value, we observe that the gen¬ 
eration of the small-scale field is quenched by the large-scale 
field—^producing a weak anticorrelation between small- and 
large-scale fields. This reciprocal correlation between the 
small-scale field and the global dynamo cycle is in agree¬ 
ment with observation s of the small-scale magnetic field at 
the solar surface ( e.g., iHagenaar et al.l 120031: |Jin et al.l l20 1 It 
I.Tin & Wa ngll2012 .12015al iFaurobert & Ricorjl2015h . In par- 
ticular, |jin et al.l ( 201 ll) have analyzed separately the small- 
scale elements based on their magnetic flux contents and the 
number of magnetic elements and found that their cyclic vari¬ 
ations for magnetic elements with fluxes in the range ~ (3- 
30) X 10^® Mx and ~ (4-40) x 10^® Mx show anticorrelation 
and c orrelation w i th sun spots, respectively. Further analy¬ 
sis bv LTin & Wangl (120121) has revealed that this anticorrelated 
component shows a solar cycle variation different from the 
sunspot variation in time and latitude, whereas the correlated 
component shows similarities with the sunspot butterfly dia¬ 
gram. Based on our simulation results in Regime IV, and the 
observational results, we conclude that the small-scale mag¬ 
netic field observed at the solar surface (particularly, the mag- 
netic elements with fluxes < 3 x 10^® Mx in the analysis of 
l.lin et al.ll201 Ih are a consequence of both the small-scale dy¬ 
namo generated locally in the solar convection zone and the 
shredding and tangling of the large-scale magnetic field gen¬ 
erated by the global dynamo. On the other hand, the magnetic 
field with a flux above 4 x 10 ^^ Mx, but smal ler than the active 
region flux in the analysis of LTin et akldlonh . are possibly the 
result of shredding and tangling of the large-scale magnetic 
field generated by th e glo bal dynamo only. Our simulations 
in Regime II (Section lT^ . when the dynamo number is large 
but the Reynolds number is small, illustrate this situation be¬ 
cause in this case the local dynamo is not possible and the 
small-scale field is still generated by the tangling of the large- 
scale magnetic field. 

^ In the upper layer of the Sun, which is our region of interest, we have 
iim ~ 10®—much larger than what we consider. However, based on stellar 
observations, we may imagine that in the Sun, D is not much larger than the 
critical dynamo number nee ded to excite the large-scale dynamo; see Sec¬ 
tion 5 of IKarak et all <2015tifl . 
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We reiterate that we observe a weak anti-correlation be¬ 
tween small-scale field and large-scale magnetic cycle only 
when the large-scale field is much stronger than the equipar- 
tition value. Based on our knowledge, we can say that the 
large-scale field in the solar convection zone may be larger 
than the equipartition value, although probably not by much. 
Therefore we expect that the small-scale magnetic field in the 
Sun is largely unaffected by the global large-scale cycle, but 
not completely. The lack or weakness of a variation of the 
small-scale magnetic field observed in our simple turbulent 
dynamo simulations certainly support the idea that both the 
small-scale dynamo and the tangling of large-scale field are 
responsible for the observed small-scale magnetic field in the 
Sun. In fact, in certain regimes our results are in agreement 
with observations of the small-scale magnetic field. However, 
future studies using more realistic solar dynamo models will 
be useful for more quantitative comparisons. 
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